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Stellar dynamics of blue compact galaxies* 
I. Decoupled star-gas kinematics in ESO400-G43 
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Abstract. In this letter we present the first results from a program aimed at measuring the stellar kinematics of 
blue compact galaxies by observing the near-infrared Calcium triplet. We show the first results for ESO400-G43 
based on deep VLT/FORS2 spectroscopy. The instabilities found in the central gaseous velocity field are not seen 
in the solid body stellar rotation curve, indicating that stars and gas are kinematically decoupled in this galaxy. 
Even if this galaxy has a perturbed gaseous velocity field, the stellar velocity dispersion on average agrees well 
with that derived from the nebular lines. 
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1. Introduction 

There has been strong interest in Blue Compact Galaxies 
(BCGs) ever since the pioneering study of Searle & 
Sargent (1972). The hypothesis that BCGs represent gen- 
uinely young galaxies, presently forming their first gener- 
ation of stars, has now been disproved for all but a few 
extreme cases (see, e.g., Kunth & Ostlin 2000 for a re- 
view). However, the evolutionary history of BCGs, their 
connection to other kinds of galaxies, and the mechanism 
that triggers the active star formation are still elusive. 

Most studies of the nature of BCGs have been based on 
surface photometry (e.g. Papaderos et al. |T996) or spec- 
troscopic analysis of nebular abundances (e.g. Izotov and 
Thuan 1999), and relatively little is known about their 
dynamics, though Hi velocity fields have been obtained 
for some galaxies (e.g. van Zee et al. 1998). The stellar 
component has been much less analysed than the gaseous 
one because the spectra are dominated by hot stars and 
ionised gas emission, diluting stellar absorption features in 



the optical domain and making kinematical and chemical 
analysis difficult. 
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* Based on observations collected at the European Southern 
Observatory, Paranal, Chile, under observing programme 71. B- 
0602. 



In Ostlin et al. <|1999l I2()()l|l we investigated the veloc- 
ity fields of the ionised gas in BCGs, using the CIGALE 
Fabry-Perot interferometer (Amram et al. 11991(1 and tar- 
geting the Ha line. These results show that many luminous 
BCGs (LBCGs) have peculiar Ha velocity fields and ir- 
regular morphology, suggesting that their starbursts have 
been triggered by dwarf galaxy mergers. Many LBCGs 
also rotate too slowly to support the stellar mass inferred 
from deep optical and near-IR photometry (Ostlin et al. 
2001). This mass discrepancy can be resolved if it is as- 
sumed that the width of the Ha line traces virial motions 
(see Melnick et al. I1987fl . indicating that these systems 
may be supported by velocity dispersion, i.e. random mo- 
tions. This is rather unexpected for such gas- rich galax- 
ies, but in line with the merger interpretation. Moreover, 
a few galaxies show peculiar rotation curves that initially 
rise rapidly, followed by a Keplerian (or even faster) de- 
cline, strengthening the interpretation that these galaxies 
are not in dynamical equilibrium (Ostlin et al. 2001). 



It is, however, well known that the ISM may be subject 
to other types of motions than purely gravitational. Winds 
and outflows, for example, are expected to be present in 
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actively star forming galaxies like these. Possibly, the Ha 
velocity field traces feedback rather than the gravitational 
potential. Even if gravitational motions are likely to dom- 
inate bulk motions on the scale of a whole galaxy, this 
possibility underlines the need for information on the kine- 
matical properties of the stellar component. 

The near-PR [Can] AA8498, 8542, 8662 triplet (hereafter Ca 
triplet) is well-suited for dynamical studies in galaxies 
dominated by young stellar populations (Dressier [1984). 
Even in a pure starburst, red supergiants will in most 
places contribute to an observable Ca triplet strength, 
and its use in probing complex galaxy kinematics has 
been demonstrated in many papers (e.g. Prada et al.[l996 
Garcia-Lorenzo et al. 119971 Kormendy and Bender 119991 
Marquez et al. I2()03jl . We have therefore initiated a pro- 
gramme aimed at studying the stellar kinematics of BCGs, 
in particular those with peculiar Ha kinematics, by using 
the Ca triplet. Here, we are interested in studying the stel- 
lar kinematics as a function of radius, meaning that levels 
fainter than the night sky background need to be reached. 
This is a challenging task due to the many sky lines and 
strong fringing of most CCDs in this wavelength region. 
However, the red-optimised FORS2 spectrograph at the 
ESO VLT is ideally suited for this work. 

Kobulnicky and Gebhardt l|1997[) studied the inte- 
grated dynamics of late type galaxies, among them a few 
BCGs, by looking at the [Oii]a3727 and Can H and K lines. 
They found that stars and gas tend to have similar global 
profiles, but their results for low mass galaxies were incon- 
clusive. Moreover, no LBCG was included in their study. 
The Ca H and K lines are not useful in our case due to 
strong emission lines in the blue spectral region. 

Galaxies with masses and properties similar to local 
LBCGs are common at higher redshift and make a signif- 
icant contribution to the cosmic star formation rate (e.g. 
Guzman et al. 1997, 2003). Hence, local LBCGs may give 
insights into the early evolution of galaxies, and have lu- 
minosities high enough to be observable even at high red- 
shift. However, the study of absorption lines relies on a 
strong continuum, and at high redshifts it will be chal- 
lenging to observe anything but emission lines. As these 
are sensitive to dynamical feedback from star formation, 
it is an important task to investigate how well-coupled 
gaseous and stellar kinematics are in local galaxies of the 
same kind. 

In this letter we report on the first results of our study: 
The first determination of the stellar kinematics of a BCG 
through the Ca triplet. In future papers we shall make a 
more detailed analysis and present the results for other 
galaxies in our sample. 

1.1. ESO400-G43 

ESO400-G43 is a luminous (M B = -19.6, H = 75 
km s _1 Mpc -1 ) BCG. It was once proposed as a young 
galaxy candidate (Bergvall & ,Torsater ll988(l , but deep sur- 
face photometry has revealed a rather massive underlying 
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Fig. 1. FORS2 I-band image of ESO400-G43 with the 
slit position indicated. The inset in the upper right corner 
shows the Ha velocity field from Ostlin et al. (1999) 

population (Bergvall & Ostlin 2002). It has a companion 
galaxy at a projected distance of 70 kpc, too far for tidal 
forces to explain the strong starburst seen in ESO 400-G43 
(Ostlin et al. lSUuT)) . 

The Ha velocity field is (upper right inset in Fig. 0, 
at first sight rather regular, but after an initial rise to 
60 km s _1 at a radius of 1 kpc, the rotation curve falls 
off rapidly, in contrast to the photometric mass profile 
(Ostlin et al. 2001). The fall is in fact faster than the 
Keplerian case, suggesting dynamical disequilibrium. The 
apparent rotational energy falls short of the photometric 
stellar mass by an order of magnitude. At much larger 
radii, ~ 15 kpc, the Hi velocity is again close to the central 
60 km s _1 (Bergvall & Jorsater Il988fl suggesting that a 
dynamical instability is responsible for the rapid decline 
of the optical rotation curve. 

We identify three possible explanations for the ap- 
parent discrepancy between dynamical and photometrical 
mass. Gas-star decoupling may be important, i.e. a stellar 
velocity field different from the gaseous one, arising per- 
haps from outflows triggered by supernova winds or infall 
in a merging process. Alternatively, global gravitational 
instabilities might mean that neither stellar nor gaseous 
motions trace the potential. Finally, the low rotational 
velocities (but not the shape of the rotation curve) could 
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Fig. 2. Two-dimensional spectrum of ESO400-G43. The 
Paschen emission lines have been marked and labelled 
with their number in the Paschen series above the spec- 
trum. The Can triplet lines and an Oi emission line have 
been marked below the spectrum. On the right, the refer- 
ence row is indicated. The height of the displayed part of 
the spectrum is 24". 
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Fig. 3. An example of subtraction of Paschen lines. In this 
spectrum, the full drawn line shows the subtracted spec- 
trum, which has also been interpolated across the emission 
lines of [Fe n] and [CI n] The dotted line shows the un- 
subtracted spectrum. 

be explained if velocity dispersion dominates the gravita- 
tional support. 

The Ha velocity field has a well-defined kinematical 
centre and position angle, and was therefore suitable for 
a study with a single slit position. Our objective was to 
derive the stellar rotation curve and velocity dispersion 
out to a radius > 2 kpc, making it necessary to reach 
surface brightness levels of /ij > 21 mag/arcsec 2 . 

2. Observations and reductions 

The observations of ESO 400-G43 were obtained between 
June and August 2003 using the spectrograph FORS2 at 
the VLT. We observed one single slit position with a posi- 
tion angle of 45° (see Fig. QJ, chosen to coincide with the 
kinematical major axis derived from the Ha velocity field 
(see Fig. □ and Ostlin et al. 1999). We used the 1028z 
grism (covering A = 773 to 948 nm) and a slit width of 
0.7", resulting in a spectral resolution of R ~ 4000 near 



the Ca triplet. The spectra were acquired as a series of 
600s exposures. In between each exposure, the telescope 
was offset along the slit. Using the same setup, six late- 
type reference stars were observed, to serve as templates 
for the cross correlation. 

We used a somewhat non-standard reduction tech- 
nique, in order to ensure effective subtraction of the strong 
atmospheric background emission at these wavelengths: 
Adjacent pairs of exposures were subtracted from each 
other, each time scaling the subtracted frame so that resid- 
uals from the OH and O2 lines were minimised. Only 
those frames resulting a in a good sky subtraction were 
retained for production of the final spectrum. Wavelength 
calibration was carried out in two dimensions using the 
OH sky lines on the unsubtracted frames. Each pairwise- 
subtracted frame was divided by a dome flatfield. The 
two-dimensional spectra were shifted to a common origin 
and coadded. Since the offset step along the slit was differ- 
ent in each case, negative residuals from objects inciden- 
tally falling in the slit cancelled when averaging the final 
frame. The resulting two dimensional spectrum (Figure[5J 
was based on 14 frames with a total exposure time of 8400 
seconds, which was used for the subsequent analysis. From 
this final two-dimensional spectrum, one-dimensional ex- 
tractions were made, and flux-calibrated using spectra of 
the standard stars LTT 7379 and LTT 7987. 

Our aim was to cross-correlate the extracted spectra 
with reference stars of different types, mainly late- type 
giants and supergiants. The galaxy spectra contain also 
narrow emission lines, most of which belong to the Hi 
Paschen series. For each of the Ca triplet lines, there is a 
Paschen line (line number 13, 15 and 16 in the series) with 
only 3-4.5 A higher rest wavelength. As the Ca absorption 
features have intrinsic widths of this order, the Paschen 
lines had to be removed before the Ca lines could be used 
to probe the stellar velocities. 

To achieve this, we first measured the fluxes and ve- 
locities of the strongest unblended Paschen lines from our 
spectra, and compared with the predicted strengths for 
Case B recombination at a temperature of 10 4 K and den- 
sity 100 cm~ 3 (Osterbrock 1989, Hummer & Storev HU87|l . 
Within the observational uncertainties, the lines follow the 
predicted values for zero reddening, consistent with the 
low extinction reported by Bergvall & Ostlin (2002). We 
then modelled the Paschen series from Pa 9 to Pa 20 for 
each extracted spectrum, scaled to match the strongest 
clean Paschen lines and subtracted the modelled emission 
spectrum (Figure|2J). The redshifts and widths of the mod- 
elled lines were determined from the strongest Paschen 
lines in the spectrum. 

Finally, we corrected for remaining strong emission 
lines like O 1 A8446, [CI 11] A8579 and [Fe 11] A8617, bad 
pixels and areas of poor sky subtraction by linearly inter- 
polating across them. The spectra, now free of emission 
lines, were cross-correlated with the reference stars, using 
the iraf task fxcor (see Tonry & Davis E2HI>- In ad- 
dition to the Ca triplet, these spectra also contain other 
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Fig. 4. Variation of heliocentric velocity (bottom) and ve- 
locity dispersion (top) in ESO 400-G43, measured from the 
absorption line spectrum. Squares with error bars repre- 
sent stellar absorption data. Crosses (x) represent PalO- 
lines, and circles [Sin] lines. Plus signs (+) show the Ha 
velocities for the same slit angle extracted from Fabry- 
Pcrot data (Ostlin et al. 1999) using a sector opening an- 
gle of 5°. The scale in kpc (bottom) and arcsec (top) is 
indicated; the zero point coincides with the reference row 
in Fig. 



weak stellar absorption features which we also included in 
the cross correlation. 

We calibrated the determination of , the line of sight 
stellar velocity dispersion, by broadening the reference 
star spectra with Gaussians of various widths, then mea- 
suring the resulting FWHM of the cross-correlation peak. 

3. Results and discussion 

In Figure 01 we show (filled squares) the stellar velocities 
and velocity dispersions along the slit, and the correspond- 
ing data for ionised gas from the PalO line (crosses) and 
[Siii]a9069 (circles). The gas velocity dispersions have been 
corrected for the instrumental width. In addition the Ha 
velocities extracted from the Fabry-Perot interferometric 
data (Ostlin et al. 11999(1 are shown (plus signs). For clarity, 
no error bars are plotted for the gas, but uncertainties arc 
in general much smaller than for the stellar component. 

The velocities for the ionised gas agree well with the 
previous investigations by Bergvall and Jorsater (Q.988), 



and Ostlin et al. (1999). The differences between the Pa 
lines and the Ha can be attributed to the poorer spatial 
resolution of the latter data. The velocities rise to ±55 
km s _1 (uncorrected for inclination) over a linear scale of 
±1.2 kpc. Outside these radii the velocities decline faster 
than the Keplerian case. On the approaching side the gas 
velocities level out at radii > 2kpc. However, the difference 
with respect to the maximum velocity is still less than 
the Keplerian prediction. Hence, the ionised gas is not in 
dynamical equilibrium. 

The stellar velocities behave differently. Firstly, the ve- 
locity gradient is much flatter than for the ionised gas. At 
r = —2.2 kpc there is a drop in the rotational velocity 
(Fig. UJ, but at r = —2.5 kpc the rotation speed is again 
high. These two outer points are though the most uncer- 
tain ones, and there is no strong indication of a decline 
similar to what is seen in the gas. On the approaching 
side, the outer points suggest a modest velocity decrease 
at r > 2.5 kpc (Fig. 0J), where a twist is seen in the Ha 
velocity field (see Fig.^ an d Ostlin et al. H999f) as if a spi- 
ral arm were present. Taken together, the stellar velocity 
field is consistent with solid body rotation over the central 
~12" (5 kpc) of our spectrum, where the S/N is good. 

The difference between the minimum and maximum 
stellar velocity is 124 km s _1 , but this is biased by the 
uncertain point at r = — 7". A straight line fit to the stellar 
velocities produces a velocity of ~ 80 km s _1 on a scale 
of 4.5 kpc (12"). With an assumed inclination of i = 55° 
(Ostlin et al. I1999|) this implies a rotational amplitude of 
50 km s _1 . Using the same simple formula as in Eq. 1 of 
Ostlin et al. IJT333I see also Lequeux 1983) with / = 0.8 
gives a mass of 1.3 x 10 9 M Q , somewhat lower than the 
photometric mass inside this radius of x 10 9 Mq. The 
photometric mass profile indicates that the rotation curve 
should continue to rise out to a radius of 4 kpc (see Fig. 5 
in Ostlin et al. I2001f) . At r = 2 kpc the predicted circular 
speed for rotational support is 801^5 km s^ 1 , whereas we 
measure 40/sin(t) km s _1 . 

This rotational measure assumes that the kinemati- 
cal position angle for the stars coincides with that for 
the gas (45°). However, the major-axis position angle 
from the broad band isophotes and Hi data is closer to 
20°. Moreover, the broad-band data suggest a smaller in- 
clination than 55°, while the Hi kinematics imply i sa 
50° (Bergvall and Jorsater IT5SS|) . Ha/H/3 mapping (Ostlin 
& Bergvall 0,994) show that the southeastern side of the 
disk is least affected by extinction and thus is probably 
nearest to us. With this orientation of the disk, an out- 
flow would have the effect, which we also observe, that 
the kinematical position angle of the gas would be ro- 
tated anti-clockwise with respect to the stellar component. 
Consequently the mass from the stellar rotation curve may 
be underestimated, and should be increased by a factor of 
1.9 (for i = 55°), or 2.3 if the true inclination is as low 
as 40°. In both cases, the mass inferred from the stellar 
rotation would be in fair agreement with the photometric 
mass. 
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The velocity dispersion measured from the ionised gas 
agrees well with that derived for the stars. Both are how- 
ever lower than the 49 km s _1 derived from Fabry-Perot 
intcrferometry (Ostlin et al. 2001). The widths we present 
in Fig. 4 are, however, close to the instrumental width and 
may be underestimated. The agreement between stars and 
gas suggests that the widths of the hydrogen recombina- 
tion lines reflect motions in response to the gravitational 
potential (see Melnick et al. 1987). Under this assumption, 
Ostlin et al. (2001) derived a mass of 8 x 1O 9 M . Taking 
cr* w 40 km s _1 we estimate the mass supported by stellar 
velocity dispersion to be 5 x 10 9 M Q . 

It is clear that the gaseous velocity field does not reflect 
dynamical equilibrium. A model that can reproduce the 
main features in both the stellar rotation curve and the 
Ha velocity field reasonably well is one with solid body 
rotation with a position angle of 20°and an outflow with 
a 30 km s _1 expansion velocity out to a radius of ~ 3 kpc, 
where the expansion is halted in a shock. However, this 
model is not well constrained. A two-dimensional stellar 
velocity field would be required in order to construct a 
physical model for the kinematics and distinguish between 
the possible explanations given in Sect. 1.1. 



4. Conclusions 

The gaseous rotation curve of ESO400-G43 increases 
steeply within 1 kpc from the kinematical centre, then 
drops faster than the Keplerian prediction. This is not 
seen in the stellar rotation curve which presents a solid 
body shape out to r = 2 kpc. Hence, stars and gas are 
dynamically decoupled in this galaxy. The reasons for this 
could be outflows due to feedback from star formation or 
dynamical instabilities in a merging process, and will be 
discussed in future papers. 

The stellar and gaseous velocity dispersions agree, 
but the amount of stellar velocity dispersion is some- 
what smaller than implied by the photometric mass of 
12±\ 6 x 1O 9 M (Ostlin et al. 2001). Nevertheless, even al- 
lowing for uncertainties in the inclination, the stellar ve- 
locity dispersion makes a larger contribution than rotation 
to the gravitational support. This is consistent with the 
picture that the starburst host galaxy could be of, or will 
evolve into, elliptical type (see Bergvall & Ostlin 2002). 

Even though the stellar and gaseous kinematics ap- 
pear decoupled in this galaxy, the mass inferred from the 
line widths of nebular lines and stellar absorption features 
agree. It remains to be seen if this is a general feature 
of luminous BCGs, but if so, it would lend credibility to 
mass estimates of high redshift galaxies that are based on 
the widths of narrow nebular lines. 
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